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Abstract. We present extensive new spectroscopy and imaging of PN G 135. 9-f 55.9. We use these data as con- 
straints to photoionization models to derive limits on the oxygen abundance. We find that PN G 135.9-1-55.9 has 
an oxygen abundance less than 1/50 of the solar value. Our models favour a value of 12-|-log O/H between 5.8 and 
6.5 de x, confirming that PN G 135.9-1-55.9 is the most oxygen-poor planetary nebula known (Tovmassian et al. 
200 1[ ). We also derive Ne/O = 0.5±0.3, S/O < 0.094, and Ar/O < 0.23. Ahhough the value of Ne/O is nominally 
high, it need not imply that the progenitor of PN G 135.9-1-55.9 converted any of its initial oxygen abundance to 
neon. The helium abundance appears to be very low, He/H ~ 0.08, but a precise determination will require a 
much more detailed study. We find that Ha/H/3 is lower than expected and perhaps variable, a finding for which 
we have no clear explanation. 

Key words, planetary nebulae: individual: PN G 135.9-1-55.9 



1. Introduction 

Recently, SBS 1150+599A has been recognized as a plan- 
etary nebula in the Galactic halo by Tovmassian et al. 
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( poo ID and renamed PN G 135.9+55.9. The spectra then 
available for this object were quite unusual for a plane- 
tary nebula, presenting only the Balmer lines of hydrogen. 
He II AA4686,5411, and very weak [O III] A5007 (~ 4% of 
H/3). A photoionization model analysis showed that such 
a spectrum implies a strongly density bounded and ex- 
tremely oxygen-poor nebula ionized by a very hot star. 
The oxygen abundance was estimated to be less than 1/50 
of the solar value, and probably between 1/100 and 1/500 
of solar assuming canonical properties for the central star, 
making of PNG 135.9-1-55.9 by far the most oxygen-poor 
planetary nebula known, with an oxygen abundance sim- 
ila r to t he lowest mea sured to date in stars (Boesgaard et 
al |1999|; Howard et al |l997[). 



In this paper, we report detailed follow-up observa- 
tions, aimed at providing more stringent constraints on 
the nature of this exceptional object. Section 2 presents 
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Table 1. Log of the spectroscopic observations 





CT)1\ ^ 1 




Ljt Hi 


WHi 


sFM3 


brMi 


date 


22 Jan 2001 


23-24 Jan 2001 


3-4 Mar 2001 


11 May 2001 


5 Mar 2002 


8 Apr 2002 


CCD 


Thomson 


Thomson 


EEV 


Tektronix 


SITe 


SITe 


format 


2048 X 2048 


2048 X 2048 


2048 X 4500 


1024 X 1024 


1024 X 1024 


1024 X 1024 


pixel size'' 


Ufj,m, o'.'ei 


14 ^m, O'.'ei 


13.5 Aim, 0'.'283 


24 /im, 0'.'36 


24 Aim, 1'.'05 


24 Aim, 1'.'05 


gain (e"/ADU) 


0.5 


0.5 


1.8 


1.4 


1.3 


1.3 


read noise (e~) 


4.8 


4.8 


3.1 


4.6 


8 


8 


spectrograph 


B&C 


B&C 


MOS 


ISIS 


B&C 


B&C 


grating'' 


300/4550A 


600/4550A 


B400/5186A 


R158R/6500A 


300/4550A 


400/5150A 


slit width 


3'.'8 


3'.' 8 


5" 


1" 


2'.'3 


2'.' 3 


spectral resolution"^ 


12. 7A 


5.6A 


23.0A 


7.8A 


8.lA 


5.8A 


wavelength interval 


3700-6700A 


4745-7200A 


3400-8000A 


6820-9740A 


3600-7100A 


3660-6760A 


arc lamp 


HeAr 


HeAr 


HgNeAr 


CuNe+CuAr 


HeAr 


HeAr 


standard stars 


HD93521 


G191B2B 


Feige 66 


Hz44 


BD+33°2642 


HD93521 




BD+33°2642 


HD93521 








BD-K33°2642 


total exp. time 


4500 s 


7800 s 


7200 s 


3600 s 


5400 s 


12600 s 


number of spectra 


3 


4 


4 


3 


3 


7 



^ These observations were obtained through clouds. 

^ Both the physical pixel size and the angle subtended on the sky are given. 

The gratings are described by their ruling (lines/mm) and effective blaze wavelength. For the grism used at CFHT, the 
ruling and zero deviation wavelength are given. All of the observations were obtained in first order. 

This is the spectral resolution measured at Ha, defined as the FWHM of the Hq line, except for the WHT spectrum, where 
this is the FWHM of P9A9229. 



the new spectroscopic data, while Section 3 deals with 
narrow-band imaging. In Section 4, we present an updated 
photoionization model analysis, taking full advantage of 
the constraints provided by our new observational data. 
This leads to a limit on the oxygen abundance which is 
now independent of any assumption about the evolution- 
ary status of the central star. In Section 5, we estimate 
the abundances of the other elements. Section 6 presents 
a brief concluding discussion. 



2. Spectroscopy 

2.1. Observations 

Table ^ presents a summary of our new spectroscopic 
observations of PN G 135.9+55.9. This table includes 
the dates of the observations, the instrumental config- 
uration, and the flux and wavelength standards that 
were used. The only observing run that suffered from 
non-photometric conditions was that of 5 Mar 2002, 
when significant cloud cover affected observations of both 
PNG 135.9+55.9 and the standard star. 

The spectroscopy from the Observatorio Astronomico 
Nacional in San Pedro Martir, Baja California, Mexico 
(SPM) was obtained using the Boiler & Chivens spectro- 
graph (B&C) and three different gratings during four ob- 
serving runs. For the 2001 observations, a rather wide slit 
(3'.'8) was used to better measure the total fluxes, while, 
for the 2002 observations, a narrower slit was used to ob- 
tain higher spectral resolution and better sensitivity to 



fainter lines. The standard stars were observed with an 
even wider slit (9"). In all cases, the slit was oriented east- 
west on the sky. Spectra of the illuminated dome wall were 
obtained to serve as flat field images. Bias images were ob- 
tained at the beginning and end of the night. 

The spectroscopy at the Canada-France-Hawaii 
Telescope (CFHT) was obtained with the Multi-Object 
Spectrograph (MOS; Le Fevre et al. [1994D . Both the ob- 
ject and the standard star were observed through a 5" slit. 
However, the observations of PN G 135.9+55.9 were ob- 
tained the night before those of the standard star. Spectra 
of the internal halogen lamp were obtained to serve as flat 
field images. 

The spectroscopy at the William Herschel Telescope 
(WHT) was obtained using the red arm of the ISIS spec- 
trograph. The object was observed through a 1" slit, while 
the standard stars were observed with a 10" slit. For these 
observations, the slit was oriented at the parallactic angle. 
Spectra of the internal lamp were obtained to serve as flat 
field images while spectra of the sky were used to correct 
for the slit illumination. 

All of the spectroscopy was reduced using the Image 
Reduction and Analysis Facility (IRAF)[| software package 
(specifically the specred package). In all cases, the over- 
scan bias was subtracted from each image. For the SPM 



' IRAF is distributed by the National Optical Astronomical 
Observatories, which is operated by the Associated Universities 
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Table 2. Raw line intensities relative to H/3 for the optical spectra'' 
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A SPMl SPM2 CFHT SPM3'' SPM4 



[Nev 


3426 






< 9.6 






[On] 


3727 






< 2.1 






[Nelll] 


3869 






1.04 ±0.52 






Hi 


3889 






2.92 ±0.73 






Hi 


3970 


10.0 ± 2.5 




6.11 ±0.64 


4.9 ±2.9 


4.5 ± 2.1 


Hi 


4101 


21.4 ± 2.8 




20.6 ± 1.2 


19.3 ±3.7 


21.3 ± 2.8 


Cii 


4267 






< 1.1 




< 3.0 


Hi 


4340 


41.0 ± 1.6 




42.05 ± 0.98 


43.8 ± 2.4 


40.7 ± 1.4 


[Oiii] 


4363 


< 1.9 




< 0.88 




< 1.8 


Hen 


4686 


76.1 ±2.3 




78.6 ± 1.5 


73.0 ±2.1 


77.8 ±1.3 


[Ar iv]'^ 


4711 


< 2.0 




< 1.0 




< 1.4 


Hi 


4861 


100.0 ± 2.1 


100.00 ± 0.84 


100.0 ± 1.6 


100.0 ± 1.7 


100.00 ± 0.85 


[Oiii] 


5007 


2.7 ± 1.3 


2.17 ±0.61 


2.87 ±0.85 


3.8 ± 1.1 


2.31 ±0.41 


Hen 


5412 


5.68 ±0.72 


5.47 ±0.62 


5.23 ±0.34 


5.64 ±0.91 


5.53 ± 0.46 


He I 


5876 


< 0.47 


< 0.69 


< 0.28 


< 1.4 


< 0.62 


Hi 


6563 


314.7 ±7.6 


296.8 ± 2.4 


260.3 ± 4.4 


288.0 ±4.1 


276.6 ± 2.0 


[Nil] 


6583 


< 3.2 


< 1.4 


< 0.45 


< 1.3 


< 0.74 


[Sii] 


6716 








< 1.3 


< 0.74 


[Sii] 


6731 








< 1.3 


< 0.74 






1.82 ±0.03 


1.33 ±0.01 


2.55 ±0.03 


0.99 ±0.01 


1.14 ±0.01 






68.4 ± 1.4 


61.4 ±0.7 


65.3 ±1.0 


61.5 ±1.3 


53.5 ±0.5 


^^(Hq)" 




646 ± 58 


508 ± 24 


541 ± 32 


494 ± 10 


514 ± 10 


Hq range' 




304-318 


256-350 


258-267 


265-324 


241-321 



When no uncertainty is given, the value represents a 2a upper limit to the flux in the line. 
^ These observations were obtained through clouds. 
These upper limits also apply to [Arlv]A4740. 

This is the total flux in emission measured at H/3 in units of 10"^* ergs~^ cm~^. 
These are the equivalent widths in A of H/3 and Ha in emission. 
' This is the total range spanned by the values of Hq/H/3 among the individual spectra. 



data, the overscan-subtracted bias images obtained dur- 
ing the night were combined and subtracted from all of the 
images. Next, the pixel-to-pixel variations were removed 
by division of the flat field image. For the WHT data, 
the slit illumination correction was then applied. The sky 
emission was subtracted during the extraction of the one 
dimensional spectra by defining sky regions on both sides 
of the object spectra and interpolating between them with 
a straight line. The wavelength calibration was performed 
using arc lamp spectra obtained at the time of the object 
observations. Finally, the spectra were calibrated in flux 
using the observations of the standard stars (Table |l|) to 
determine the instrumental sensitivity function. The in- 
dividual spectra were calibrated in both wavelength and 
flux before being summed together. 

Table |^ presents the raw line intensities relative to 
H/? measured in the optical spectral region for PN G 
135.9±55.9. The line intensities presented in Table ^ are 
those for the summed spectra from each observing run. 
The line intensities were measured using the software de- 



Table 3. Raw line fluxes for the WHT spectra 



scribed by McCall et al. (1985). The uncertainties quoted 
for each line intensity (Icr) include contributions from the 
fit to the line itself, from the fit to the reference line, and 



ion 


A 


flux'' 


F(A)/f(H/3)'' 


[Ar III] 


7135 


< 0.45 


< 0.06 


Hi 


8750 


5.7 ±2.3 


0.74 ±0.31 


[Sin] 


9069 


< 0.59 


< 0.08 


Hi 


9229 


20.3 ±3.4 


2.64 ±0.45 



The fluxes are given in units of 10"^^ ergs"^ cm~^. When 
no uncertainty is given, the value is a 2a upper limit to the 
flux in the line. 

These flux ratios are relative to F(H/3) measured for the 
CFHT spectrum on a scale where /(H/3) = 100. 



from the noise in the continuum for both the line and ref- 
erence line. When only a limit is given, this corresponds 
to a 2(7 upper limit to undetected lines. 

Table ^ presents the raw fluxes and the intensity ra- 
tios relative to H/3 for the WHT spectrum. The only line 
deflnitely detected is P9A9229; P12A8750 is detected at 
only the 2a level. Again, when no uncertainty is given. 
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Fig. 1. We compare the spectra of PNG 135.9+55.9 and the standard star G191B2B on an arbitrary magnitude scale. 
For PNG 135.9+55.9, we plot the CFHT spectrum for A < 7500A and the WHT spectrum for A > 6760A, neither 
corrected for reddening. The WHT spectrum wa s nor malized to the CFHT spectrum as de scribe d in the text. For 
G191B2B {B-V = -0.32 mag), we plot the Oke (\l99(\) fluxes. As noted by Tovmassian et al (|200l[) , PN G 135.9+55.9 
has a remarkably blue continuum. 



the value represents a 2a upper limit to the line inten- 
sity. These line intensities and limits were measured using 
IRAF's splot routine. The fluxes represent the fluxes mea- 
sured directly in the summed WHT spectrum. The inten- 
sity ratios relative to H/3 were computed adopting the H/3 
flux from the GFHT spectrum and correcting the WHT 
fluxes for the difference in the slit widths used. Based upon 
the spatial profile of Ha from the CFHT spectrum, a 5" 
slit intercepts 3.32 times more nebular emission than a 1" 
slit. The WHT fluxes were then multiplied by this fac- 
tor when computing the relative intensities presented in 
Table ^. Comparing the continuum fluxes measured in the 
CFHT and WHT spectra in the 6760-7500A region, where 
the fringing in the CFHT spectrum is not too severe, the 
continuum flux in the WHT spectrum should be scaled 
upwards by a factor of 2.07 to match that in the CFHT 
spectrum. This scale factor is in good agreement with the 
value of 1.95 expected based upon the 1" slit used for 
the WHT spectra and the 1'.'22 seeing measured from the 
spatial profile of the continuum in the summed spectrum. 



Generally, there is excellent agreement among the line 
intensities over the wavelength range H(5-He iiA5412. In 
the near- ultraviolet, the SPM spectrograph has very low 
efficiency and the upper limits we derive from those spec- 
tra are considerably less restrictive than the detections or 
limits from the CFHT spectrum. We give upper limits to 
the [S ii] line intensities only for the last two SPM runs, 
since the CCD used for the first two SPA^I observing runs 
and that at CFHT suffered from fringing in the red. 

The notable exception to the good agreement among 
the line intensities is Ha. There is significant variation 
in the Ha/H/3 ratio among the summed spectra for the 
different observing runs and between individual spectra 
for at least the SPM2 and SPM4 observing runs. In both 
the SPM2 and SPM4 data sets, the dispersion among the 
Ha fluxes for the individual spectra also significantly ex- 
ceeds that for the H/3 fluxes. In the last line of Table ^ 
we indicate the range of Ha/H/3 values found among the 
individual spectra during each observing run. This vari- 
ation is very puzzling, since we normally obtained all of 
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Table 4. Log of the direct imaging observations 
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Telescope 


Date 


CCD 


Instrument 


Fiher'' 


Exposure time'' 


SPM 2.1m 
NOT 2.6m 


27-28 May 2001 
1-2 June 2001 


SITeS'^ 
CCD7'* 


Mexman 
ALFOSC 


Ha (6565A, llA) 
red cont. (6650A, 46A) 

y#18 (5470A, 220A) 
Ha #21 (6564A, 33A) 


3900s (5) 
1500s (3) 

2700s (4) 
2700s (4) 



The central wavelength and the bandpass width for each filter are given in parentheses. 
^ The number of images is given in parentheses. 

This CCD has 24 /im pixels in a 1024 x 1024 format. Its gain and readnoise are 1.3e~/pix and 8e~, respectively. The plate 
scale is O'.'312/pix. 

This CCD has 15 /im pixels in a 2048 x 2048 format. Its gain and readnoise are ^ le~/pix and 6e~, respectively. The plate 
scale is 0'.'l88/pix. 



the spectra consecutively on the same night (SPM2 is the 
exception). If this variation is real, it is occurring (irreg- 
ularly) on a time scale of the order of an hour. Such be- 
haviour is not at all expected in a nebular plasma (e.g., 
AUer 1987). In a typical SPM spectrum (of a half hour 



duration), at least 50,000, 15,000, and 2,500 photons are 
detected at Ha, H/3, and H7, respectively, so the vari- 
ation we see in the Ha/H/3 ratio would not appear to 
be due to poor photon statistics. Over this wavelength 
range, we do not see any variation exceeding more than 
a few percent in any of the standard star observations. 
We made no effort to orient the slit at the parallactic 
angle, but the wide slits used, particularly for the 2001 
observing runs, should compensate for the effects of dif- 
ferential refraction. Regardless, were differential refraction 
the culprit, we should see nearly equally larg e vari ations 
in H7/H/3 as we see in Ha/H/3 (Fihppenko |l982| ), but 
we see none. Three instrumental effects, however, affect 
the 2001 data from SPM. First, the object was acquired 
by blind offset, so the centering of the object in the slit 
was almost certainly not optimal. Second, the offset guider 
is known to flex relative to the instrument field of view, 
so the object centering was likely somewhat variable for 
the 2001 observations at SPM. Finally, the spectrograph 
was out of focus due the CCD being mis-aligned with the 
camera's focal plane. It is not clear, however, how any 
of these might introduce variations in the Ha/H/3 ratio 
alone without affecting other line ratios. None of these is- 
sues affect the CFHT data nor those from SPM in 2002, 
yet the Ha/H/3 variations exist in these data sets as well. 
Tovmassian et al. (2001) found similar variations, from a 
variety of observing sites, though they attributed them to 
the poor observing conditions affecting their observations. 
Although unusual, it would appear that the variations in 
Ha/H/3 are real. 

In Fig. ^ we compare the spectra of PNG 135.9+55.9 



and the standard star G191B2B (Oke |1990|). For 
PNG 135.9+55.9, we plot the CFHT spectrum for A < 
75OOA and the WHT spectrum for A > 6760A, without 
applying any reddening correction to either spectrum. The 



WHT spectrum was scaled upwards by a factor of 2.07, 
as described previously. As noted by Tovmassian et al. 
(2001), this planetary nebula has a remarkably blue con- 
tinuum. 

Finally, these new data do not provide any direct di- 
agnostic of the physical conditions in the nebular plasma. 
No density diagnostic has been detected to date, though 
[Ar iv]AA471 1,4740 might have been expected given the 
high degree of ionization. Similarly, the only temper- 
ature diagnostic available is the upper limit to the 
[O III] A4363/5007 ratio and it provides no useful constraint 
unless the density is unusually high, 10^ — 10^ cm~'^, which 
is excluded given the nebular flux and size (see Sec. 4). 



2.2. Reddening 



From the Schlegel et al. (199S) reddening maps, the ex- 
pected foreground reddening is E{B — V) ^ 0.01 — 0.02 
mag. The H7/H/3 and H5/H/3 ratios imply E{B - V) 
0.3 — 0.35 mag, based upon a temperature of 2 x 10"^ K, a 



density of 10 — 10 cm ^, the Storey fc Hum mer (1995) 



line emissivities, and the Fitzpatrick (1999) monochro- 
matic reddening law parametrized with a ratio of tota l- 
to-sclective extinction of 3.041 (McCall & Armour ^OOOl ). 
On the other hand, the Ha/H/3 ratio implies a reddening 
E{B ~ V) < 0.23 mag for the same physical conditions, 
reddening law, and line emissivities, even if we consider the 
largest line ratio we observe, Ha/H/3 = 3.5. For Ha/H/3 
ratios at the low end of the range observed, the reddening 
is zero or negative. The intensity of P9/H/3 from the WHT 
spectrum implies E{B — V)= 0.05 mag. We can also com- 
pute a reddening using He iiAA4686,5412. Adopting the 
same physical conditions, reddening law, and line emis- 
sivities, we find negative reddenings, i.e., HeiiA4686 is 
too bright relative to HeiiA5412 by 12%. In any case, 
it appears that the reddening is at most modest, with 
E{B — V) < 0.3 — 0.35 mag. In the remainder of this pa- 
per, we shall assume that the reddening is zero. None of 
the conclusions would be affected had we adopted a mod- 
est reddening. 
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3. Imaging 

3.1. Observations 

We obtained direct images of PN G 135.9+55.9 in narrow- 
band filters from two sites (Table ^). Images were obtained 
in the emission line of Ha and a nearby continuum band- 
pass with the 2.1m telescope in SPM. The Mexman filter 
wheel was used for observations. The seeing was modest, 
with stellar images having FWHM of 1'.'5 or worse. Images 
were also obtained using ALFOSC at the 2.6 m Nordic 
Optical telescope (NOT) on La Palma, Canary Islands, 
Spain. A narrow filter centred at Ha was used, but the 
continuum was measured using a Stromgren y filter. The 
NOT images benefited from better seeing, « 1'.'2, on the 
first night of this run. The spectrophotometric standard 
stars HZ44 and Feige 67 were observed along with the ob- 
ject on the second night for calibration purposes. During 
both imaging runs, sequences of at least 3 images in each 
filter were obtained, with one long exposure usually not 
exceeding 20 min and two shorter ones. This tactic allowed 
us to keep cosmic ray rate within reasonable limits in or- 
der to later combine the images for improved statistics 
and to permit the removal of cosmic rays. 

The images were processed using IRAF. The steps in- 
cluded bias substraction, flat field correction, image com- 
bination, and apperture photometry. The STSDAS pack- 
age within IRAF was used for the deconvolution of the 
images. 

The primary uses of the Ha images were to determine 
the angular size of the nebula and to derive the density 
distribution of the gas. Two selected NOT images from the 
first night, where the seeing was better and the focus was 
good were combined to produce the image presented in Fig 
^. The nebula appears to be slightly elliptical in shape. No 
substructure is evident in the main body of the nebula. 
Its size comfortably exceeds the PSF of stellar objects, 
so the instrumental profile should not affect its shape. 
Nonetheless, we conducted series of deconvolution exper- 
iments on the nebul ar images usi ng th e Lucy-Richardson 
method (Lucy 1974 Richardson 1972 ) that showed that 
the size and shape of the nebula were unchanged after de- 
convolution. This is our strongest evidence that the cen- 
tral star does not emit significantly in Ha. For comparison 
the image of the central star obtained in the Stromgren y 
filter is presented in the bottom panel of Fig. ^ with the 
Ha contours of the nebula superposed. We used the SPM 
observations, where the continuum filter was much closer 
in wavelength to the line filter, to substract the image of 
central stellar object from that of the nebula. Figure || 
presents the spatial profiles of the nebular emission in Ha 
from the NOT image without subtracting the central star 
and from the SPM image after substraction of central star. 
Except for a slight depression around the peak, the shape 
of the nebula in both observations is very similar. 

The NOT observations in the Stromgren y filter were 
used to estimate the magnitude of the central star. We 
inferred — 17.9 ± 0.15 mag, corresponding to a flux 



02 
40'00" h 
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Fig. 2. The upper panel presents the image of 
PN G 135.9+55.9 in the light of Ha, while the lower panel 
presents the Stromgren y image of the central star with 
the Ha contours of the nebula superposed. The continuum 
and Ha images were aligned using field stars. Both images 
were obtained with the NOT. 



of ^5556 = 2.52 X 10~^^ ergcm"^ s~^ A in a good agree- 
ment with the flux determined from spectroscopy. This 
continuum measurement should be reliable since there are 
almost no lines in the spectrum of this planetary nebula, 
the contribution of HeiiA5413A being negligible. 

3.2. The nebular angular diameter and density 
distribution 

The angular diameter of the planetary nebula is an impor- 
tant parameter since, in combination with the H/3 flux, it 
allows a determination of the electron density and the neb- 
ular mass. We followed the recipe of Bedding & Zijlstra 



(1994) to assure that our measurements correspond to 
common standards. (For cautionary remarks, see van Hoof 
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Angular Size (arc sec) Normalized Angular Size (t^^nax^'^-^ '^^^ sec) 



Fig. 3. This figure presents the spatial profiles of the neb- 
ula from the Ha images and of the central object from the 
6650A continuum image. The filled and open squares rep- 
resent cuts along major and minor axes of the nebula from 
the SPM images after subtraction of the continuum. Open 
triangles represent the spatial profile along the major axis 
of the Ha image obtained at the NOT without continuum 
subtraction. The dotted line is the profile of the central 
star (SPM). The fluxes are normalized to the maximum 
intensity of the SPM continuum-subtracted image, except 
for the central star, which is normalized to its maximum 
intensity. 



2000| .) There are a variety of definitions of the diameter 
of a PN in the literature: the FWHM, the surface brigth- 
ness contour at 10% of the maximum surface brightness, 
and the outermost contour at which emission is found. 
Since PN G 135.9-1-55.9 appears to have a slightly elliptical 
shape in the Ha images, we have measured the diameter 
of the nebula along two axes. We find that the diameter of 
the nebula is ll'.'4x9'.'8 if we use the outermost contour, 
6'.'3x5'.'2 if we use the 10% surface brightness contour, and 
3'.'5x2'.'95 if we use the FWHM. 

Figure ^ shows the density and flux distributions 
adopted in our models. The symbols correspond to the 
observed Ha emission profiles across the major and mi- 
nor axes obtained from the SPM images (shown in Fig. 
^). The curve corresponds to a spherical model with a 
density distribution described by n = exp — (0/2.8)^, 
where 9 is the angular distance to the center in seconds 
of arc. In this model, the emissivity of Ha is taken to be 
simply proportional to ti^ and is computed for a uniform 
temperature. Such a simple representation accounts very 
satisfactorily for the observed Ha emission profile. 



Fig. 4. The density distribution used for the model calcu- 
lations is presented in the lower panel. The predicted flux 
distribution from our models is drawn with a solid line in 
the upper panel, supposing that the emissivity is propor- 
tional to and computed at a constant temperature. The 
open and filled squares denote the spatial distributions of 
the nebular emission along the major and minor axes, re- 
spectively, obtained from the SPM observations, as shown 
in Fig H The abscissa denotes relative angular size nor- 
malized to the maximum, azimuthally-averaged radius of 
5'.'6. 



In the analysis of the photoionization models that fol- 
lows, we adopt a nebular radius of 5'.'6, i.e., similar to the 
diameter of the outermost contour. 



4. Photoionization analysis: new limits on the 
oxygen abundance 

4.1. The modelling procedure 

The observational data do not permit a direct estimate 
of the oxygen abundance, since nothing is known of the 
electron temperature or of the presence of oxygen ions 
more charged than 0++. It is therefore necessary to rely 
upon photoionization models. 

We have constructed sequences of photoionization 
models in which the oxygen abundance varies over sev- 
eral orders of magnitude. The models are constrained by 
the available observations, which consist of the observed 
line intensities, the equivalent width of H/3, the total fiux 
in H/3, the size of the nebula, and the radial distribution 
of Ha shown in Fig. ^. As already noted by Tovmassian et 
al. (2001), the shape of the stellar continuum only implies 
that the star is hotter than 50,000 K. The models are com- 
puted with the photoionization code PHOTO, using the 
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Table 5. Some global properties of the sequences of models investigated 



Q(H0) L{HI3) d z Rout M„e6 

[K] [erg s"^] [erg cm"^ s"^] [kpc] [pc] [cm] [ Mq] [cm"^] 



100 000 


6 10-*^ 


0.48 


25.0 


20.7 


1.8 10^** 


0.78 


125 


100 000 


1 10*" 


0.080 


10.2 


8.5 


6.3 10" 


0.070 


190 


125 000 


6 10'*" 


0.29 


19.2 


15.9 


1.6 10^® 


0.45 


148 


125 000 


1 lO"*" 


0.048 


7.8 


6.5 


6.1 10" 


0.044 


223 


150 000 


6 10*" 


0.20 


16.0 


13.3 


1.3 10^** 


0.29 


163 


150 000 


1 10*" 


0.032 


6.5 


5.4 


4.8 10" 


0.027 


245 



atomic data listed in Stasinska & Leitherer (1996). The 
central star is assumed to radiate as a blackbody of tem- 
perature Ti,. The hydrogen density at a radius r is taken 
to be n = ric cxp —{r/hY', where ric is a free parameter 
and h = 2.8d/2.06 10^, where d, the distance to the nebula 
(in the same units as h and r), is also a free parameter. 
The ionizing radiation field is treated in the outward-only 
approximation. The computations start close to the star 
and are stopped when the equivalent width in H/3, W(H/3), 
becomes equal to the observed value, taken to be 70A. 

With such a representation of the nebula, it is easy to 
show that, for each assumed stellar temperature T^, and 
total luminosity L^,, the total observed flux in H/3 implies 
a certain distance to the planetary nebula. The observed 
angular size of the image then fixes the outer radius of 
the nebula. By trial and error, we determine the value of 
He for which the radius Rout corresponding to W(H/3) = 
70A is close to 5.6d/2M 10^. Table | gives some global 
properties of selected sequences of photoionization mod- 
els. This includes the distance, d, and the height above 
the Galactic plane, z, implied by these models. Note that 
Rout and the nebular mass, Mneb, are decreasing func- 
tions of the electron temperature. The values given in the 
table correspond to the low metallicity end of the model 
sequences. 

The chemical abundances in each sequence of mod- 
els are parametrized by 0/H, with the abundances of 
the heavy elements with respect to oxygen following the 
recipe of McGaugh (1991). Note that the abundance ra- 
tios in PN G 135. 9-f 55.9, as in any PN, may be signifi- 
cantly different from those assumed in McGaugh's recipe, 
especially for carbon and nitrogen. For helium, we as- 
sume an abundance of 0.08 with respect to hydrogen in 
all models. This is close to the value estimated from the 
observed Hell 4686/H/3 ratio and using the upper limit 
for Hei A5876 given by the CFHT spectrum (see Sect. 5). 
In any event, our estimate of the chemical composition 
of PN G 135.9+55.9 is independent of the relative abun- 
dances adopted in the model sequences since our estimate 
uses the observed line intensities and relies on the ion- 
ization and temperature structure of the nebula, which 
depend essentially on hydrogen and helium in the domain 
of interest. 



Fig. |5| presents the results of our models as a func- 
tion of 12 + log(0/H). Each row of panels corresponds 
to a different central star temperature: = 100, 000 K, 
r^= 125, 000 K, and = 150, 000 K, as indicated in the 
first panel. Each column of panels displays a different line 
ratio. In each panel, two series of models are represented 
with different symbols. Circles correspond to models with 
central stars having a total number of ionizing photons, 
Q(II°), equal to lO"*^ photons sec~^, while squares corre- 
spond to models with a total number of ionizing pho- 
tons of 6 X lO**^ photons sec~^, values that roughly bracket 
the luminosities of the central stars of planetary nebulae 
as computed by post-AGB evolutionary models (Blocker 
1995). The horizontal lines show the observed values: thick 
lines for measured values or limits, thin lines indicating 
the uncertainties in measured values, and upward- and 
downward-pointing arrows denoting lower and upper lim- 
its, respectively. For the observational data, we adopted 
the CFHT observations, since they provide the most ac- 
curate measurements and the most stringent limits. For 
the Ha/H/3 ratio, however, we plot the value derived from 
the SPMl observations. As noted earlier, our observations 
give inconsistent results for the Ha/H/3 ratio. Values of 
Ha/H/3 around 3.1 are easily accounted for by our mod- 
els, but lower values are extremely difficult to explain. We 
will return to this issue in Sect. 4.3. 

4.2. Discussion of tlie models 

The new observations allow us to eliminate the models 
with > 125, OOOK, since these produce a [Ne v]/[Ne III] 
ratio much larger than is observed. Our models show that 

should be around 100,000 K. T^, cannot be much lower 
than this value since He i A5876 would then be observed. 
Note that the observed lower limit to [O iii]/[0 ii] pro- 
vides no useful constraint upon the ionization structure 
(or central star temperature) of our models. 

For the models with T-^ = 100, OOOK shown in Fig. 
|, the oxygen abundances compatible with the observed 
[O iii]/H/? ratio are 5.8 < 12+logO/H< 6.3 dcx. Allowing 
for a reasonable uncertainty in implies 12 + log 0/H < 
6.5 dex. The upper limit to [Ne v]/[Ne ill] eliminates mod- 
els excited by a star with a temperature significantly above 



Richer, M. G. et al.: The a-element abundances in PN G 135.9+55.9 



9 




log (0/H) 



12 + log (0/H) 



12 + log (0/H) 



12 + log (0/H) 



12 + log (0/H) 



5 6 7 
12 + log (0/H) 



Fig. 5. This figure presents our grid of photoionization models. Each row of panels corresponds to a particular central 
star temperature, indicated in the first panel. Each column of panels plots a different line intensity ratio. In each panel, 
two sequences of models are included: circles denote models with ionizing luminosities of 10^^ photons s~^, squares 
denote models with ionizing luminosities of 6 x 10'*^ photons s^^. The thick horizontal lines indicate the observed values 
or their limits from the CFHT spectrum. In the case of limits, upward- and downward-pointing arrows denote whether 
they are lower or upper limits, respectively. Thin horizontal lines denote the uncertainty range of measured values. 




5 6 7 8 9 
12 + log (0/H 



5 6 7 
12 + log (0/H) 



Fig. 6. Here, we compare two models with — 150, 000 K and Q(II'^) = lO^^photons^^, but with the circles denoting 
a filling factor of 1 (the model from Fig. |^) and squares denoting a filling factor of 0.1. As this experiment demonstrates, 
the maximum oxygen abundance allowed by the models decreases as the filling factor decreases. 



100, 000 K which would otherwise permit 12 + log 0/H of 
the order of 7-8 dex. 

We have constructed other series of models to test the 
robustness of these conclusions. For example, we have cal- 
culated models in which the nebular radius, total nebular 



flux, and H/3 equivalent width were varied, but the con- 
clusions remain unchanged. 

We have also considered models in which the gas is dis- 
tributed in small clumps with the same global density law 
as above, but with an overall filling factor of 0.1. Although 
the Ha image is smooth at our resolution, we cannot a pri- 
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ori exclude the presence of small scale clumps or filaments. 
Clumpy models that fit the observational constraints will 
result in a lower global ionization level than smoothly- 
distributed models. For the purposes of illustration, Fig. || 
compares the smooth model of Fig. |with = 150, 000 K 
and (5(H") = lO^'^phsec"^ (circles) with a clumpy model 
whose filling factor is 0.1 (squares). Because the ioniza- 
tion is lower in the clumpy model, the 0+++/0++ ratio 
is lowered with respect to the smooth model, while the 
[O III] A5007/H/3 ratio is raised. Consequently, the oxygen 
abundance compatible with the observed [O ill] A5007/H/3 
ratio is smaller than in the smooth case. Note that the 
model presented here gives 12 -I- logO/H ~ 6.5 dex, but 
still violates the [Ne v] A3426/[Ne iii] A3869 and Ha/H/3 
constraints. The total nebular mass of a model of given 
total H/3 luminosity and given radius is roughly propor- 
tional to the square of the filling factor. Therefore, clumpy 
models that satisfy the observational constraints will have 
lower nebular masses than the corresponding models listed 
in Table ^. Finally, if the density in the clumps were ex- 
tremely high, models with high metallicities, even as high 
as solar, could account for the weak intensities of the for- 
bidden lines, because these lines would be quenched by 
collisional de-excitation. For this to occur, densities ex- 
ceeding 10® cm~'^ are required for both [O ill] A5007 and 
[Ne III] A3869. Given that the H/3 flux and equivalent 
width as well as the size of the nebula are known, such 
high densities would imply that the filling factor would 
have to be of the order of 10~^ or less, which is highly 
unrealistic. 

Our computations were made assuming that the star 
radiates as a blackbody. However, a more realistic stellar 
atmosphere would give a different spectral energy distri- 
bution for the ionizing photons, particularly at the largest 
energies. One expects that extended, metal-poor atmo- 
spheres could have larger fluxes at energies above 100 eV. 
This would increase the [Ne v]/[Ne ill] ratio and conse- 
quently strenghten our conclusion that 12 -I- logO/H < 
6.5 dex. On the other hand, absorption by metals could 
depress the number of photons able to produce [Ne v] 
(Ranch 2002), in which case a star with - 150,000K 
could become acceptable. In this case, however, the ex- 
citation of the nebula would be lower than that pre- 
dicted by blackbody models with Ti,= 150, 000 K and 
the line ratios would resemble those produced by the 
blackbody model with = 100, 000 K, again implying 
12 -I- logO/H < 6.5 dex. In any case, the amount of met- 
als in the atmosphere is not expected to be large, unless 
the atmosphere contains dredged-up carbon. A definitive 
answer to this problem can only come from a direct mea- 
surement of lines from more highly charged ions. 

Note that the distance we obtain for our object (see 
Table ||) indicates that it is located in the Galactic halo, in 
agreement with its radial velocity (Tovmassian et al 2001). 
Its derived nebular mass is compatible with the range of 
nebular masses derived for Magellanic Cloud PNe (Barlow 
1987). 



To conclude this section, we emphasize that our 
new observational data allow us to confirm that 
PN G 135.9-1-55.9 is an extremely oxygen-poor planetary 
nebula, with an oxygen abundance less than 1/50 of the 
solar value. Our models favour a value of 12 -I- logO/H 
between 5.8 and 6.5 dex, compared with the solar value of 
12 -f logO/H = 8.83 dex (Grevesse & Sauval p^98| ). Our 
modelling experiments indicate that this conclusion is in- 
dependent of plausible changes in the properties assumed 
for the central star or nebular envelope. 



4.3. The Ha/Hf) problem 

Our preferred models with T^, = 100, 000 K and low metal- 
licities are compatible with the Ha/H/3 ratio derived from 
the SPMl data, as well as with the highest values found 
in individual spectra from the other SPM observations 
(Table ||, Fig. ^). However, as discussed earlier, we find ap- 
parently significant variations of the Hq;/H/3 ratio between 
the different observing runs and even among individual 
spectra obtained during individual runs. The majority of 
the individual spectra indicate that Ha /H/3 is belo w 3, a s 
do the observations reported by Tovmassian et al (2001). 

Because of collisional excitation of the hydrogen lines, 
Ha/H/3 cannot have the recombination value, but is ex- 
pected to be larger. Collisional excitation is unavoidable 
at electron temperatures above ~ 15, 000 K as soon as 
there is a small fraction of residual neutral hydrogen. The 
amount by which Ha/H/3 exceeds the recombination value 
(2.75 at 20,000 K, 2.70 at 30,000 K, using the case B coef- 
ficients of Storey & Hummer 1995) depends upon both the 
electron temperature, which is higher for higher values of 
Ti,, and on the amount of neutral hydrogen. For example, 
the models with — 125, OOOK predict a value of Ha/H/3 
around 3.4 at low metallicities. With our observational 
constraints, there is not much room for a drastic reduc- 
tion of the proportion of neutral hydrogen in our models. 
Indeed, for a given and chemical composition, the pro- 
portion of H'^ at each point in the nebula is completely de- 
termined by (and roughly proportional to) ne/{Lv/R^), 
where ng is the local electron density, Ly is the stellar 
luminosity in the V band and R is the distance of this 
point to the star. Ly jB? is a distance-independent quan- 
tity that relates the stellar flux in the V band to the an- 
gular distance of this point to the star; is determined 
by the density distribution law obtained from our Ha im- 
ages and the value of ric is imposed by the observational 
constraints on W(H/3) and the total angular radius, as ex- 
plained in Sect. 4.1. Our models of PN G 135.9-^55.9 with 

= 100, OOOK indicate that the Ha/H/3 ratio should not 
be lower than 3 if the nebula is metal-poor. Values of 
T-t significantly below 100,000 K could be consistent with 
some of the observed values of Ha/H/3 but they are ex- 
cluded by the failure to observe He i A5876 in this object. 

Docs this mean that the object is not as oxygen- 
poor as inferred above? There are several arguments 
against a higher abundance. In our T^, = 100, OOOK mod- 
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els, Ha/H/3< 3 implies 12 + logO/H > 7.5 dex, but 
this is clearly incompatible with the observed value of 
[O III] A5007/H/3, as seen in Fig. |. At the other ex- 
treme, for our T^, = 150, 000 K models, Ha/H/3 < 3 implies 
12 + logO/H > 8 dex, which, though marginally compat- 
ible with the observed [O ill] A5007/H/3, predicts a value 
for the [Ne v] A3426/[Ne ill] A3869 ratio far larger than 
observed. 

Nor is it likely that any possible variability of the 
Ha/H/3 ratio affects our oxygen abundances. Supposing 
that this variation is real and refers entirely to the nebu- 
lar radiation, one would then expect the [O ill] A5007/H/3 
ratio to vary as a consequence of variable ionization or 
temperature conditions, but this is not seen. All of the line 
intensities apart from Ha remain remarkably constant, in- 
cluding the [O III] A5007 line. Therefore, we believe that 
the Ha/H/3 problem does not affect our conclusions as 
regards the oxygen abundance in PNG 135.9+55.9. 

However, this Ha/H/3 problem is extremely puzzling. 
One does not expect this ratio to vary in nebular con- 
ditions for an extended object. One possibility could be 
that PN G 135.9+55.9 harbours a compact disk. Such 
a suggestion has been made for other planetary nebu- 
lae based upon either morphological, spectroscopic, or 
variability arguments (He 2-25: Corradi 1995; IC 4997: 
Miranda & Torrelles 2000, Lee & Hyung 2000; and M 2-9: 
Livio & Soker 2001). Accretion disks in close, interact- 
ing binary systems have the particularity of both being 
variable (Warner 1995| ) and having Ha/H/3 ratios much 
smaller than 3, sometimes attaining values below unity 
(e.g., Williams 1995). If PN G 135.9+55.9 contained such 
an accretion disk and if this disk contributed to the emis- 
sion of the hydrogen and helium lines in the central part 
of the nebula, this could explain both the variability of 
the Ha/H/3 ratio in PN G 135.9+55.9 and the values of 
2.7 or lower observed in some of our spectra. However, an 
accretion disk would be an unresolved point source in our 
images, and our deconvolution experiments found no sig- 
nificant contribution to the Ha emission from the central 
source. Likewise, the lack of Dopper-broadened emission 



erence. The second column of panels gives the value of 
r(0++) defined as: 



lines (Tovmassian et al 2001) also argues against an ac- 
cretion disk as the or igin o f a significant fraction of the 
line emission (Warner 1995 ). 



5. The abundances of the remaining elements 

Our high signal-to-noise spectra enabled us to measure 
the intensity of the [Ne III] A3869 line, and to derive upper 
fimits to the intensities of [S ill] A9069 and [Ar ill] A7135, 
which allow us to make some inferences on the abundances 
of Ne, S and Ar. 

To compute abundances, we must rely upon the elec- 
tron temperature and ionization structure provided by our 
models. The relevant quantities are shown in Fig. ^for the 
models from Fig. ^. The first column of panels on the left 
repeats the [O iii] A5007/H/3 ratio from Fig. || for ref- 



r(o 



/ Ten{0++)nedV 
J n{0++)nedV 



(1) 



The next columns of panels give the ionization frac- 



tions x(0++)/a;(H+), a;(0++)/x(Ne++), x{0++)/x{S-^ 
and a;(0++)/a;(Ar++), where the notation a;(0++) stands 
for: 



x(0^ 



/ n{0++)nedV 
Jn{0)nedV ' 



(2) 



The value of a;(0++)/a;(H+) is shown here simply to 
demonstrate that O^^ is far from being the dominant 
oxygen ion in the nebula, and that a more robust determi- 
nation of the oxygen abundance would require observing 
lines from higher ionization stages. 

Using the same atomic data as in the photoionization 
code and taking as a characteristic temperature of the 
emission of the [Ne ill] A3869 and [O ill] A5007 lines a 
value of 30,000 K derived from our preferred models (Fig. 
I), we find that Ne++/0++= 0.70 ±0.4, where the uncer- 
tainty reflects the uncertainties of the [O ill] A5007 and 
[Ne III] A3869 line intensities. We furthermore note that 
both 0++ and Ne++ are minority ions (Fig. |^). Generally, 
one adopts Ne+'^/0++= Nc/0. In our preferred models 
(those corresponding to Ti, = 100,000 K), we find from Fig. 
Rthat the ionization correction factor is rather around 0.7, 
thus leading to Ne/0 = 0.5 ± 0.3. 

In a similar fashion, our observed limits on the intensi- 
ties of the [Ariv]AA4711,4740 lines allow us to derive that 
Ar+++/0++ < 0.45. From our models (Fig. 0) we find 
that the ionization correction factor is 0.5, and therefore 
Ar/0< 0.23. 

Given that the nebula is strongly density-bounded, the 
observed limits on the [N ii] A6584 and [S ii] A6716, A6731 
line intensities give only crude limits on abundance ratios 
involving these elements. With the ionization correction 
factors from our preferred model, we find N/H < 2.3x 10~^ 
and S/H < 7x 10^^. Such limits are not very useful, except 
to infer that the N/H ratio is at most equal to that in the 
Orion nebula. 

Our infrared spectra allow us to estimate upper limits 
to the intensities of [S ill] A9069 and [Ar ill] A7135, which 
are respectively 8 x 10"'* and 6 x 10~* of the intensity 
of H/3. These upper limits imply that S++/0++ < 0.031 
and Ar++/0++ < 0.012. In our preferred photoionization 
models, the ionization correction factor to derive Ar/0 
from Ar^^/0++ is around 20, so that the upper limit 
on [Ar III] A7135 implies Ar/0 < 0.24 (in agreement 
with the upper limit derived from [Ar iv]AA471 1,4740). 
Regarding sulfur, we find that the ionization correction 
factor is around 3 for S/0 from our models (Fig. im- 
plying that S/0 < 0.094. It must be noted that the atomic 
data concerning the ionization structure of S and Ar are 
uncertain (see e.g. Ferland et al. 1998). Our models, were 
computed without including dielectronic recombination to 
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Fig. 7. This figure presents various physical parameters for the grid of models in Fig. ||. From left to right, the columns 
present [O ill] A5007/H/3, the temperature in the 0^"+ zone, the weighted fraction of oxygen in the form of and 
the ionization correction factors for Ne++, S"'"^, and Ar+++ with respect to O'*'"''. These data are required to derive 
the abundances of Ne, S, and Ar in PN G 135.9+55.9. 



low- lying levels for these ions. It is likely that the real ion- 
ization fractions of S++ and Ar++"'" are actually higher 
than predicted by the models, giving smaller ionization 
correction factors and more stringent limits. 

It is interesting to compare the ratios of Ne/0, S/0 
and Ar/0 we find for PN G 135.9+55.9 with those of other 
kinds of objects. Table || shows the values for a sample of 
PNe in the Galactic disk (Kingsburgh & Barlow 1994) and 
in the Galactic halo (Howard et al. 1997), for the Orion 
nebula (Esteban et al. 1998), and for the Sun (Grevesse 
& Sauval 1998). This table is of course subject to uncer- 
tainties. Even the solar abundances are quite uncertain 
for Ne, Ar and O (see Grevesse & Sauval 1998; the recent 
oxygen abundance determination from AUende Prieto et 
al 2001 yields a value that is only 73% that obtained by 



Table 6. a-element abundance ratio comparison with 
other objects 



Grevesse & Sauval 1998). One might then argue that the 
Ne/0 ratio in PNG 135.9+55.9 is compatible with the so- 
lar value, but we note that it is about twice the value 
found in the Orion nebula and in disk planetary nebu- 
lae, where the systematic errors in the abundance deriva- 
tions are likely to be similar. Supposing our Ne/0 ratio 
is corect, it might indicate some conversion of O into Ne 



object 


Ne/0 


S/0 


Ar/O 


PNG 135.9+55.9 


0.5 + 0.3 


< 0.094 


< 0.23 


Galactic disk PNe 


0.26 


0.017 


0.005 


Galactic halo PNe " 


0.13 


0.016 


0.0016 


Orion nebula 


0.18 


0.03 


0.014 


Sun " 


0.18 


0.03 


0.004 



This work. 

Kingsburgh & Barlow (1994). 

Howard et al. (1997). These authors also show that the 
abundance ratios are more dispersed in halo PNe than in 
disk PNe. 

Esteban et al. (1998). 
Grevesse & Sauval (1998). 



by a capture. A few similar cases are known among plan- 
etary nebulae (e.g., BB-1 has a Ne/0 = 0.77; Howard et 
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al. 1997). Another possibility is that the material from 
which the progenitor of PNG 135.9+55.9 formed had an 
anomalous Ne/0 ratio. Abundance studies of metal-poor 
stars indicate that the very early galaxy was chemically- 
inhomogeneous, with individual sites of star formation 
being influenced by the explosions of nearby supernovae 
(e.g., Burris et al. 2000). The yields of O and Ne from 
individual supernovae are also a function of the stellar 
mass (Woosley & Weaver |1995| ; Thielemann et al |1996D , 
while observations of metal-poor halo stars indicate that 
the scatter in oxygen abundances is 0.3-0.5 dex at very 
low oxygen abundances (e.g., Israelian et al 2001). As a 



result, it is probably not surprising that the progenitor of 
PN G 135.9+55.9 might not have formed out of material 
with the same Ne/0 ratio that has characterized the more 
recent Galactic disk. Regardless, any conversion of O into 
Ne is at most modest. Even were all of the Ne formed by 
nuclear processing from oxygen, the initial oxygen abun- 
dance would have been only 0.18 dex larger than our pre- 
ferred values. Thus, the extremely low oxygen abundance 
in PN G 135.9+55.9 is genuine and not due to nuclear and 
mixing processes in the progenitor star. The oxygen abun- 
dance in PN G 135.9+55.9 should consequently reflect the 
chemical composition of the medium out of which the star 
was made. The limits we obtain on S/0 and Ar/0 in our 
object are consistent with this view (although the limits 
are not very stringent). 

The accuracy of the helium line intensities achieved 
in the present observations raised the hope of obtaining 
very accurate helium abundances. Taking the intensities 
determined from our CFHT observations and adopting a 
temperature of 30,000 K as inferred from our models and 
case B coefficients from Storey & Hummer (1995) wc find 
that He++/H+= (7.50 ± 0.14) x 10"^. The upper limit 
to the intensity of He i A5876 in the same spectra gives 
an upper limit to He+/H+ of 0.24 x 10"^ when using the 
emissivities from Benjamin et al. (1999) in the low density 
limit. A proper determination of the uncertainty in the 
derived helium abundance should, however, account for 
collisional excitation of the lines, a possible small amount 
of reddening, possible deviations from case B, possible un- 
derlying absorption, as well as temperature gradients in- 
side the nebula. All of this can only be attempted once 
the Ha/H/3 problem is solved. It is therefore premature 
to propose an accurate value for the helium abundance in 
PNG 135.9+55.9. 

6. Conclusions 

Our new, extensive observations of PNG 135.9+55.9 con- 
firm its unusual nature. Spectroscopy covering the wave- 
length interval 3400-9700A reveals lines of only Hi, Heii, 
[Oiii], and [Neiii]. Deep Ha imaging was used to derive 
the radial density distribution. 

Using these data as constraints, we constructed a new 
set of nebular models from which we derived the abun- 
dance of oxygen and the Ne/0, S/0, and Ar/0 abun- 
dance ratios. We confirm the extremely low value of the 



oxygen abundance, which we find to be less than 1/50 of 
the solar value: our models favour a value of 12 + log 0/H 
between 5.8 and 6.5 dex. The distance implied by these 
models places PNG 135.9+55.9 in the Milky Way halo, 
in accordance with its radial velocity (Tovmassian et al 
^001 ). The models also imply nebular masses in the range 
expected. 

For the a-element ratios, we find Ne/0 — 0.5 ± 0.3, 
S/0 < 0.094, and Ar/0 < 0.23. The Ne/0 ratio may 
be somewhat higher than is commonly found in plane- 



tary nebulae in the Milky Way disk (e.g., Henry 1989 



Kingsburgh & Barlow 1994). One possibility is that the 
progenitor of PN G 135.9+55.9 converted some of its O 
to Ne. It is also possible that the anomalous Ne/0 ra- 
tio is the result of discrete chemical enrichment in the 
very early evolution of the galaxy (e.g., Burris et al. 
^2000 ). Regardless, of the cause, any conversion of O to 
Ne has been modest and does not affect our conclusion 
that PN G 135.9+55.9 is the progeny of an intrinsically 
very oxygen-poor star. 

An unusual characteristic of PN G 135.9+55.9 is its 
low Ha/H/3 ratio, for which we find no clear explanation. 
Despite its low metallicity and the concomitant high elec- 
tron temperature that should result in coUisionally excited 
Balmer lines of Hi, PN G 135.9+55.9 has an Ha/H/3 ra- 
tio typically below 3. Furthermore, Ha/H/3 appears to be 
variable between observing runs and even within a single 
night. One possible explanation for both the low Ha/H/3 
ratio and its variability is if PN G 135.9+55.9 contains an 
accretion disk, though the evidence is not convincing. At 
any rate, this issue does not appear to affect our conclu- 
sions regarding the chemical abundances. 

We also measure a low He/H ratio of ~ 0.08. This 
makes PNG 135.9+55.9 interesting as a probe of the pre- 
galactic He abundance. However, the derivation of a very 
precise He abundance will require the resolution of a num- 
ber of outstanding issues, including the Ha/H/3 problem, 
the foreground reddening, and the internal temperature 
structure. 
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